Effect of proton-proton Coulomb repulsion on soft dipole excitations of light 

proton-rich nuclei 
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We perform three-body model calculations for soft dipole excitations of the proton-rich Borromean 
nucleus ^^Ne. To this end, we assume that ^^Ne takes the ^''0-fp-l-p structure, in which the two 
valence protons are excited from the O"*" ground state configuration to 1~ continuum states. We 
employ a density-dependent contact force for the nuclear part of the pairing interaction, and dis- 
cretize the continuum states with the box boundary condition. We show by explicitly including the 
Coulomb interaction between the valence protons that the Coulomb repulsion does not significantly 
alter the El strength distribution. We point out that the effect of the Coulomb repulsion in fact 
can be well simulated by renormalizing the nuclear pairing interaction. 



PACS numbers: 21.10.Gv,23.20.-g,21.60.-n, 27.20.-|-n 

Properties of unstable nuclei with large excess of pro- 
ton or neutron are one of the most important current 
topics of nuclear physics. In several radioactive beam fa- 
cilities in the world, many unstable nuclei far from the 
/3-stability line have been discovered In particular, 
neutron-rich unstable nuclei have been extensively stud- 
ied and some exotic features have been observed. This 
includes a large concentration of the dipole strength dis- 
tribution at low energies 5, 6], referred to as a "soft dipole 
excitation" . This type of excitation is naively under- 
stood as an oscillation between weakly bound valence nu- 
cleon(s) and the core nucleus The relation between 

the soft dipole excitation and a largely extended density 
distribution, that is, a halo or skin property has been dis- 
cussed for light neutron-rich nuclei, both experimentally 
[i,|i,[i-[i3 and theoretically [H-ifll- Recently, the neu- 
tron halo structure has been discussed also in the ^^Ne 
nucleus, based on the measured large Coulomb break-up 
cross sections plj[. 

In contrast to neutron-rich nuclei, proton-rich nuclei 
have been less studied. It has not been fully clarified 
whether similar exotic features are present also in proton- 
rich nuclei. For instance, a proton halo structure in e.g., 
^B, 12n, i^F, and "Ne has been discussed d [H-il], but 
no clear evidence has been obtained so far. 

In order to investigate proton-rich unstable nuclei and 
discuss their similarities and differences to neutron-rich 
nuclei, it is indispensable to assess the effect of the 
Coulomb repulsion between valence protons. In the pre- 
vious work, we analyzed the ground state properties of 
^^Ne using a three-body model [IE] ■ We have shown that 
the effect of the Coulomb repulsion is weak enough and 
the two valence protons in the ground state of ^'^Ne have 
a spatially compact configuration, that is, the diproton 
correlation, similar to a dineutron correlation in neutron- 
rich Borromean nuclei [26l - [30l ] . We have also shown that 
the effect of the Coulomb interaction between the va- 
lence protons can be well accounted for by renormaliz- 
ing the nuclear interaction. A similar conclusion was 
achieved recently also by Nakada and Yamagami, who 



performed Hartree-Fock-Bogoliubov (HFB) calculations 
for A^=20,28,50,82, and 126 isotones 

In this paper, as an continuation of the previous study, 
we discuss the effect of the Coulomb repulsion on excited 
states of ^^Ne. Our interest is to investigate whether a 
similar renormalization for the Coulomb repulsion works 
also for the soft dipole excitation, which plays an im- 
portant role in the astrophysical two-proton capture on 

We assume the ^^Ne nucleus as a three-body system 
composed of an inert spherical core nucleus ^^O and two 
valence protons. The three-body Hamiltonian in the 
three-body rest frame reads 



H = /,W+/i(2) + Pi:Pa + Fpp(ri,r2), (1) 
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+ l^pc(r»), 



(2) 



where m and Aq are the nucleon mass and the mass 
number of the core nucleus, respectively, /i^*^ is the 
single-particle (s.p.) Hamiltonian for a valence proton, 
in which /i = mAc/{Ac -I- 1) is the reduced mass and 
Vpc is the potential between the proton and the core nu- 
cleus. The third term in Eq.([T]) is a two-body part of the 
recoil kinetic energy of the core nucleus. For the proton- 
core potential Vpc, we employ a Woods-Saxon (WS) plus 
Coulomb potential in the same manner as in the previ- 
ous work |25|]. We use the same parameters for the WS 
potential as in those listed in Ref. [25[ . 

We solve the three-body Hamiltonian, Eq. ([1]), by ex- 
panding the wave function on the uncorrelated basis as, 
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where 
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(1-2) ] . (4) 

Here, (/)/cm(r) is a s.p. wave function with k = {n,l,j), 
while J and M are the total angular momentum of the 
two-proton subsystem and its z component, respectively. 
CikiM is the expansion coefficient. The summation in 
Eq. ([3]) is restricted to those combinations which satisfy 
TT = (— for a state with parity tt. 

In the actual calculations shown below, we include the 
s.p. angular momentum I up to 5. We have confirmed 
that our results do not change significantly even if we 
include up to a larger value of Z. In order to take into 
account the effect of the Pauli principle, we explicitly 
exclude the lsi/2, lp3/2: and lpi/2 states from Eq. ([3]), 
which are occupied by the protons in the core nucleus. 

For the pairing interaction Vpxii we assume a density- 
dependent contact interaction [13, E^, HI] together with 
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PP 



Vp^\ as in the 



the Coulomb repulsion, Vpj 
previous work [25]. We take a cutoff energy Ecut = 30 
MeV and include those configurations which satisfy 
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where enij is a s.p. energy [l9|. Within this truncated 
space, we determine the strength of the nuclear part of 
the pairing interaction, Vpp\ using the empirical value 
of the neutron-neutron scattering length, ami — —18.5 
fm T^. The parameters for the density dependence are 
adjusted so as to reproduce the experimental value of the 
two-proton separation energy of ^''Ne, S2p — 0.944 MeV. 

In our calculations, the continuum s.p. spectra are dis- 
cretized within a box of i?box=30 fm. Thus, energies of 
the two-proton 1~ states as well as the El strength dis- 
tribution are also discretized. The El strength function 
from the ground state is defined as 



(6) 



Using the strength function S{E), we can also calculate 
the fc-th moment of energy defined as 



Sk = / dEE''S{E) = ^(^t^OV 



(9) 



Notice that 5*0 and Si correspond to the direct and 
energy-weighted sum of dB{El)/dE, respectively. 

From the completeness of the 1^ basis, we can estimate 
the sum-rule-values as 
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where r2N-c = (ri + r2)/2 is the distance between 
the center of mass of the two valence protons and the 
core nucleus. For the core nucleus "'^^O {Aq = 15 and 
Zc = 8), we obtain S'o,sr = 1.49 e^fm^ and S'i,sr = 5.69 
e'^fm^MeV. Due to the Pauli forbidden transitions, the 
actual value of Sq is smaller than 5*0. sRj while Si is larger 
than 5*1, SR. 



TABLE I: The results for the soft dipole excitations of ^^Ne 
obtained with the three-body model of So and Si 

are the non-energy weighted sum rule and the energy weighted 
sum rule, respectively. -Bcont = Si /So is the centroid energy 
of the dipole strength distribution. ^-Bcont is a shift of the cen- 
troid energy with respect to the result of the exact treatment 
of the Coulomb interaction. 



pairing 


So 


Si 


-fi'ccnt 


SEccnt 




(e^fm^) 


(e^fm^MeV) (MeV) 


(MeV) 


Nucl. -1- Coul. 


1.206 


11.02 


9.140 





Nucl. only 


1.206 


10.45 


8.666 


-0.47 


Ren. Nucl. 


1.205 


10.86 


9.017 


-0.12 


No pairing 


1.206 


15.50 


12.86 


3.72 



Our main results are summarized in Fig. [T]and Table 
m For a plotting purpose, we smear the Dirac delta func- 
tion in the El distribution S{E) with a Cauchy-Lorentz 
function 



dB{El) ^ r 



1 



dE 



tt{e- nw,Y + r2 



(12) 



where hwi — Ei — £'g.s. , -Eg.s. being the energy of the 
ground state, and fXi is the i?(El) strength for «-th 1^ 
two-proton state. 



= 3 (f f ] ^0 ] *g.s. / , 
with the El operator given by 
'Ac-Zc 



' Ac + 2 
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[riVi^(ri) + r2Vi^(f2)]. (8) 



with the width parameter F of 1.0 MeV. In order to 
discuss the effect of the Coulomb part of the pairing 
interaction, we also perform the calculations with two 
other treatments for the Coulomb interaction. One is to 
switch off the Coulomb interaction ("Nucl. only"), that 
is, Vpp = Vpp\ keeping the same values for the parame- 
ters of the nuclear pairing interaction Vpp"^ as in the full 
calculation ("Nucl.-|-Coul."). The other is again to use 
the nuclear interaction only, but renormalize the parame- 
ters, that is, Vpp = Vpp'^ ("Ren. Nucl."). To renormalize 
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FIG. 1: (Color online) Comparison of the El strength dis- 
tributions for ^^Ne obtained with several treatments for the 
Coulomb interaction between the valence protons. The solid 
line is obtained by fully including the Coulomb interaction, 
while the dashed hne is obtained by switching off the Coulomb 
repulsion. The thin sohd line is the result of the renormal- 
ized nuclear interaction, which is readjusted to reproduce 
the ground state energy without the Coulomb interaction. 
The dotted line denotes the result without any pairing in- 
teraction. These distributions are smeared with the Cauchy- 
Lorentz function, Eq. with T = 1.0 MeV. 



the interaction, we use the empirical proton-proton scat- 
tering length, flpp ~ —7.81 fm, instead of the neutron- 
neutron scattering length ann, and determine the other 
parameters so that the two-proton separation energy S2p 
is reproduced. This leads to about 10.0% reduction of the 
strength of the pairing interaction. Notice that this value 
for the reduction factor is consistent with the finding of 
Ref. [3l|. See Ref-JHI for further details of the procedure. 
For a comparison, we also show the results without any 
pairing interaction. These treatments for the Coulomb 
interaction are applied only to the excited states, while 
the same ground state wave function is used for all the 
cases. That is, the ground state is calculated with the 
full treatment of the Coulomb interaction, that yields 
76% of ((^5/2)^ and 16% of (si/2)^. This ground state is 
used also for the no-pairing calculation for the dipole ex- 
citations. (These values for the occupation probabilities 
are slightly different from those in Ref.jl^, as we use a 
smaller i?cut in this paper. We have confirmed that the 
dipole strength distribution does not change much even 
though we use the smaller value of i?cut-) Notice that our 
results for Sq are consistent with the result of Grigorenko 
et al [13, that is. So = 1.56 e^fm^ for (51/2)^=48% and 



So = 1.07 e^fm^ for (51/2)2=5%. Table U also lists the 
centroid energy defined as i?cent = Si/ So [32,] , and its rel- 
ative value with respect to the result of the full treatment 
of the Coulomb interaction. 

From Fig. [1] and Table HI we can see that the pairing 
interaction shifts considerably the El strength distribu- 
tion towards the low energy region, similarly to the dipole 
distribution in neutron-rich nuclei |2C| . This large 
shift of the El strength distribution originates mainly 
from the nuclear part of the pairing interaction. If the 
Coulomb part is switched off, the strength distribution is 
shifted only slightly, as is shown in Fig. [T]by the dashed 
line. The shift of the centroid energy, SEccnt, is only 
— 0.47MeV. The sign and the magnitude of SEccnt is con- 
sistent with the result of Hartree-Fock-|-RPA calculations 
for medium- heavy nuclei shown in Ref. [s^, although 
SEccnt for the soft dipole excitation in ^^Ne is somewhat 
larger. 

The result with the renormalized nuclear pairing in- 
teraction is shown by the thin solid line in Fig. [T] As 
one can see, the result of the full calculation (the thick 
solid line) is well reproduced by this prescription. We 
can thus conclude that the renormalization works well 
not only for the ground state [2^, [Slj , but also for the 
dipole excitations. 

In summary, we discussed the influence of the Coulomb 
repulsion between the valence protons upon the soft 
dipole excitation in ^^Ne. We showed that the effect of 
the Coulomb repulsion is so weak that the main feature 
of the dipole response is similar between neutron-rich and 
proton-rich weakly bound nuclei. We also showed that 
the effect of the Coulomb interaction can be well mocked 
up by renormalizing the nuclear pairing interaction. This 
renormalization works both for the ground state and for 
the dipole excitations, and thus from a practical point 
of view one can use a renormalized pairing interaction in 
order to understand the structure of proton-rich nuclei. 

One of the current topics of proton-rich nuclei is two- 
proton radioactivity. Bertulani, Hussein, and Verde ar- 
gued that the final state interaction plays an important 
role in discussing the energy and the angular correlations 
in a two-proton emission process ^33;]. It would be an 
interesting future work to investigate how the renormal- 
ization works for those correlations. A work towards this 
direction is now in progress. 



This work was supported by the Grant-in-Aid for 
Scientific Research (C), Contract No. 22540262 and 
20540277 from the Japan Society for the Promotion of 
Science. 



[1] I. Tanihata et al.„ Phys. Rev. Lett. 55, 2676 (1985) ; [2] A. Ozawa et al., Phys. Lett. B 334, 18 (1994). 

Phys. Lett. B 160, 380 (1985) ; Phys. Lett. B 206, 592 [3] A. Ozawa, T. Suzuki and I. Tanihata, Nucl. Phys. A 693, 
(1988). 32 (2001). 



4 



[4] B. Jonson, Phys. Rep. 389, 1 (2004). 

[5] N. Fukuda et al, Phys. Rev. C70, 054606 (2004). 

[6] T. Nakamura et al, Phys. Rev. Lett. 96, 252502 (2006). 

[7] K. Ikeda, T. Myo, K. Kato, and H. Toki, Lecture Notes 
in Physics 818, 165 (2010). 

[8] H. Sagawa, Nucl. Phys. A 583 755 (1995). 

[9] T. Kobayashi et al, Phys. Lett. B 232, 51 (1989). 
[10] R. Arnc et al, Phys. Lett. B 250, 19 (1990). 
[11] S. Shimoura et al., Phys. Lett. B 348, 29 (1995). 
[12] T. Aumann et al., Phys. Rev. C 59, 1252 (1999). 
[13] A. Leistenschneider et al, Phys. Rev. Lett. 86, 5442 
(2001). 

[14] R. Paht et al, Phys. Rev. C 68, 034318 (2003). 

[15] P.G. Hansen and B. Jonson, Euro. Phys. Lett. 4, 409 

(1987). 

[16] C. Bertulani and G. Baur, Nucl. Phys. A 480, 615 (1988). 
[17] G.F. Bertsch and H. Esbensen, Ann. Phys. (N.Y.) 209, 
327 (1991). 

[18] H. Esbensen and G.F. Bertsch, Nucl. Phys. A 542, 310 

(1992). 

[19] H. Esbensen, G. F. Bertsch and K. Henckcn, Phys. Rev. 

C 56, 3054 (1997). 
[20] H. Esbensen, K. Hagino, P. Mueller, and H. Sagawa, 

Phys. Rev. C 76, 024302 (2007). 
[21] T. Nakamura et al, Phys. Rev. Lett. 103, 262501 (2009). 



[22] M. V. Zhukov and L J. Thompson, Phys. Rev. C 52, 

3505 (1995). 

[23] L. V. Grigorcnko, Yu. L. Parfcnova, and M. V. Zhukov, 

Phys. Rev. C 71, 051604(R) (2005). 
[24] L. V. Grigorenko, K. Langanke, N. B. Shul'gina and M. 

V. Zhukov, Phys. Lett. B641, 254 (2006). 
[25] T. Oishi, K. Hagino and H. Sagawa, Phys. Rev. C82, 

024315 (2010) ; Phys. Rev. C82, 069901(E) (2010). 
[26] K. Hagino and H. Sagawa, Phys. Rev. C72, 044321 

(2005) ; K. Hagino and H. Sagawa, Phys. Rev. C76, 

047302 (2007). 

[27] M. Matsuo, K. Mizuyama, and Y. Serizawa, Phys. Rev. 

C 71, 064326 (2005). 
[28] M. Matsuo, Phys. Rev. C 73, 044309 (2006). 
[29] N. Pillct, N. Sandulescu, and P. Schuck, Phys. Rev. C 

76, 024310 (2007). 
[30] Y. Kanada-En'yo, N. Hinohara, T. Suhara, and P. 

Schuck, Phys. Rev. C 79, 054305 (2009). 
[31] H. Nakada and M. Yamagami, Phys. Rev. C83, 

031302(R) (2011). 
[32] Tapas Sil, S. Shlomo, B. K. Agrawal and P. -G. Reinhard, 

Phys. Rev. C 73, 034316 (2006). 
[33] C. A. Bertulani, M. S. Hussein, and G. Verde, Phys. Lett. 

B666, 86 (2008). 



